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Abstract 

The stable configurations of non-rotating and rotating hybrid stars composed of colour superconducting 
quark matter core are constructed using several equations of state (EOSs). We use a set of diverse EOSs for 
the nuclear matter which represents the low density phase. The EOSs at higher densities correspond to the 
quark matter in the colour superconducting phase and are computed within the NJL-like model for different 
values of the scalar diquark and vector current couplings strengths. The phase transition to the quark matter 
is computed by a Maxwell construction. We find that the stability of the hybrid stars are mainly governed 
by the behaviour of the EOSs for the colour superconducting quark matter. However the compositions of 
hybrid star are sensitive to the EOS of the nuclear matter. The value of the critical rotation frequency for the 
hybrid star depends strongly on the EOS of the nuclear matter as well as that for the colour superconducting 
quark matter. Our results indicate that the EOS for the colour superconducting quark matter can be obtained, 
by adjusting the parameters of the NJL model, to yield the stable configurations of the hybrid star having 
the maximum mass ~ \.5Mq in the non-rotating limit and the critical rotation frequency ~ 1 kHz. 
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I. INTRODUCTION 



The present knowledge of quantum chromodynamics (QCD) suggests that quark matter might 
be in different color superconducting phases at high densities. Thus, one expects the core of the 
hybrid stars to be composed of color superconducting quark matter (CSQM) surrounded by a 
nuclear mantle. The possible CSQM phases are the two-flavor color superconductor (2SC) 1 1-3], 
the colour flavour locked (CFL) phase M |5|], and crystalline color superconductor (CCS) [Ig, lZD. 
The speculation that the CSQM exists in the core of the hybrid stars has triggered many theoretical 
investigations both on the modeling of the equation of state (EOS) of quark matter and on the 
phenomenological signatures of the presence of quark matter in the compact stars [jsl . 

The nuclear matter phase of the hybrid star is described b;^ the various models which can 
be broa dly gro uped into (i) non-relativistic potential models [9], (ii) non-re 



models lilOMlSn . (iii) field theoretica' 



Brueckner-Hartree-Fock model UJ 



based relativistic mean-field models lll4l - ll6ll and (iv) Dirac 



ativistic mean-field 



- 2011 . The quark matter in the colour superconducting phases 
are usually described either within the MIT bag model or using a more realistic NJL-like model. 
The studies based on the MIT bag model indicate the existence of stable configurations of hybrid 
stars with the CFL quark matter core [|2ll - l24ll . Further, the MIT bag model predicts the absence 
of the 2SC colour superconducting phase in the hybrid stars pSl). The scenario is some what 
different when NIL model is employed to study the hybrid stars with CSQM core. The stab 



configurations of hybrid stars with 2SC quar k matter core are possible within the NIL model [|26- 



2911 . However, earlier investigations i30l- l32|l based on the NIL model ruled out the possibility of 



CFL quark matter at the core of the hybrid stars, because, it rendered the hybrid star unstable. 



Only very recently [133 



3411 . it has been demonstrated that inclusion of the six-fermion interaction 



term together with large enough values of the scalar diquark coupling strength in the NIL model 
can yield stable configurations of the hybrid star containing 2SC or CFL quark matter core. The 



NIL model is also app 



the hybrid stars [|35 



36|1 



led to study the possibility of existence of the CCS quark matter phase in 



The stability and the structure of the non-rotating hybrid stars are quite sensitive to the choice 



of the EOS of the nuclear matter and the quark matter ||21 , 



3711 . Further, one often finds that even 



though the stable configurations of the non-rotating hybrid star for a given EOS belong to the 
third family of compact stars, but, the maximum rotation frequency upto which these hybrid stars 
are stable is much lower than the corresponding mass- shedding (Keplerian) frequency Issl . 39]. 
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The EOSs for the quark matter in the unpaired or in the various colour superconducting phases 
employed in these investigations were obtained within the MIT bag model. Recently I40I1 . a more 
realistic EOS for the unpaired quark matter computed within the NIL model is used to show that 
the maximum mass of the non-rotating hybrid stars depends sensitively on the choice of the EOS 
of the nuclear matter. It is necessary to construct stable configurations of the non-rotating and 
rotating hybrid stars using realistic EOSs for the nuclear matter and for the quark matter in the 
colour superconducting phases. 

In the present work, we compute several EOSs and use them to study the properties of the non- 
rotating and rotating hybrid stars composed of CSQM core. The lower density part of these EOSs 
correspond to the nuclear matter and are based on the variational and mean-field approaches. Our 
set of EOSs for the nuclear matter around the saturation density (po = 0.16 fm~^) is constrained by 
the bulk properties of the finite nuclei. But, their behaviour at densities, p > po, are significantly 
different. The EOSs for CSQM are calculated within the NIL model using different values for 
the scalar diquark and vector current coupling strengths. The EOS at intermediate densities are 
obtained using a Maxwell construction. 

The paper is organized as follows, in Sec. II we describe , in brief, the models employed 
to construct the EOSs for nuclear matter and the CSQM. In Sec. Ill we present the results for 
the equilibrium sequences for non-rotating and rotating hybrid stars. In Sec. IV we state our 
conclusions. 

II. EQUATIONS OF STATE 

We compute the EOSs which correspond to the nuclear matter at lower densities and CSQM in 
the 2SC or CFL phases at higher densities. The EOS at intermediate densities are obtained using 
a Maxwell construction. For nuclear matter in the jS equilibrium, we employ a set of diverse EOSs 
which are obtained using various approaches, like, variational, non-relativistic mean field (NRMF) 
and relativistic mean field (RMF). In Fig. [T]we plot various nuclear matter EOSs. The low density 
behaviour these EOSs are very much similar as they are constrained by the bulk properties of the 
finite nuclei. But, their behaviour at higher densities are so different that the resulting neutron 
star properties are at variance. In Table HI we list some key properties of the non-rotating neutron 
stars obtained using these nuclear matter EOSs. It can be seen from TableUthat the values of the 
maximum neutron star masses are in the range of 2.0 - 2.SMq and the radius Ria at the canonical 
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neutron star mass vary between 11.3 - 14.8km. It is interesting to note that the values of the 
maximum neutron star mass for both the APR and TMl EOSs are equal, but, the radius at the 
canonical mass of the neutron star is reasonably smaller for the APR EOS. This is due to the fact 
that the APR EOS is softer relative to the TMl at intermediate densities and it becomes stiffer at 
high densities as can be seen from Fig. [T] Similar is the case with SLY4 and BSRIO EOSs. We 
shall see in next section that these pairs of nuclear matter EOSs, for which the maximum neutron 
star masses are the same, yield significantly different structure for the hybrid stars. 

The EOSs for the CSQM in the 2SC or CFL phase are obtained within the NIL model. The 
input variables of the NIL model are the chemical potentials for all the quark flavours and colours 
in the chemical equilibrium which is given by the matrix: 

/ij = {^5^" + l^QQf)5ai, + [M3 (T,U + (.T,U 6"^ (1) 

where, yu is the quark chemical potential, ixq is the chemical potential of the electric charge equal 
to minus the electron chemical potential jXe and /ia and yUg are the colour chemical potentials asso- 
ciated with the two mutually commuting colour charges of the SU{3)c gauge group. The explicit 
form of the electric charge matrix Qf = diagf(^, -|, -|), and for the colour charge matrices 
Ti, = diagci^, -^,0), and VSTg = diagc{\, \, -1). In the mean-field approximation, the pressure 
at vanishing temperature reads as. 



/. = 4is:(r„(r,cr,-— ^|A,p-2G5 2cr,2, + ^ + — ^ | dkk^ \ei\ + - B (2) 

^ c=l a=l ^ i=l 

where, cr„ j are the quark-antiquark condensates and A,, are the three diquark condensates. The 
values of cr, and A<^ are determined using, 

1^=0 (3) 

OCT; 



In Eq. ^ cl>o is the mean field expectation value for isoscalar vector like meson co given as [|32h 
ojq = 2Gv{QM I il/liipu + ^li^pd + ^l^s I QM). This field modifies also the chemical potentials: 
lJ-u,d,s f^u,d,s - ojq. The 6, are the dispersion relations computed by following the Ref. [Si]. The e,- 
depend explicitly on the values of current quark masses, quark-antiquark and diquark condensates 
and various chemical potentials appearing in Eq. ([T]). The = //^/(12;r^) is the contribution to 
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the pressure from the electrons. The constant B is so determined that the pressure vanishes at zero 
density and temperature. In addition to the Eqs. ([3]) and dH), the pressure must satisfy, 

_ dp 



dfiQ 



= 0, (5) 



"3 ^ 1^ = 0, (6) 
«,.|^ = 0, (7) 

SO that local electric and colour charge neutrality conditions are met. Once, the pressure as a 
function of quark chemical potential is known, quark matter EOS can be easily computed. 

The model parameters, the current quark masses m„rf^i, quark-antiquark coupling Gs, the 
strength K of the six fermion or "t Hooft" interaction and the cutoff parameter A are taken to 
be [41], 



m„ = = 5.5 MeV, (8) 

m, = 140.7 MeV, (9) 

GsA^ = 1.835, (10) 

KA^ = 12.36, (11) 

A = 602.3 MeV. (12) 

After fixing the masses of the up and down quarks, m„ = = 5.5 MeV, the other four parameters 



are chosen to reproduce the following observables of vacuum QCD [|41ll : nij^ = 135.0 MeV, niK = 
497.7 MeV, = 957.8 MeV, and = 92.4 MeV. This parameter set gives = 514.8 MeV. 
The value of B for this set of parameters is (425.4MeV)'^. There are two more parameters, the 
diquark coupling strength Go and the vector current coupling strength Gy, which are not known. 
One expects that the diquark coupling has a similar strength as the quark-antiquark coupling. We 
construct quark matter EOS for Go = 1.1 - I.2G5 with Gy = - 0.2Gs . 

In the 2SC phase, pairing occurs only between the u and d quarks and the s quarks remain 
unpaired leading to Ai = A2 = and A3 0. On the other hand, in the CFL phase, Ai ^ 0, A2 ^ 
and A3 0. In the left panels of Figs. [2]and[3l we plot the pressure as a function of the quark 
chemical potential for the nuclear matter and for the quark matter in the 2SC and CFL phases. The 
phase realized at a given chemical potential is the one having largest pressure. Thus, it is evident 
from the P - // curves that direct transition from nuclear matter to the CFL quark matter occurs 



for the case of APR and SLY4 EOSs. For the TMl and NL3 EOSs, transition from nuclear matter 
to the CFL quark matter proceeds via 2SC phase at intermediate densities. For the BSRIO EOS, 
nuclear matter to CFL quark matter phase transition proceeds via 2SC phase only for Go = 1 -SGj 
with Gy = 0. We see from these figures that the pressure, at which the transition from the nuclear 
to quark matter occurs, decreases with increasing or decreasing Gy- For instance, pressure at 
the phase transition reduces almost by a factor of two with increase in Go from 1.1 G5 to I.2G5. 
The solid circles on the various EOSs for the nuclear matter indicate the values of Pi 4 which is the 
pressure at the center of the neutron star with the canonical mass (1.4Mq). The values of transition 
pressure is close to that of Pia for the cases plotted in the lower and the upper panels of Figs. |2] 
and[3l respectively. For the completeness, in the right panels of Figs. [2]and[3l we display the plots 
for the pressure as a function of the baryon density for the case of APR, BSRIO and NL3 EOSs. 
The phase transition to the quark matter is computed by a Maxwell construction, 
ni. HYBRID STARS WITH CSQM CORE 

We construct the equilibrium sequences of the non-rotating and rotating compact stars using the 
EOSs obtained in the last section. These EOSs correspond to the nuclear matter at lower densities 
and the CSQM in the 2SC or CFL phase at higher densities as shown in Figs. [2] and [31 The nuclear 
matter EOSs are taken from the published literature as summarized in Table HI The EOSs for the 
CSQM are computed within the NJL model for different values of the scalar diquark coupling 
strength Go and the vector current coupling strength Gy. The other parameters of the model are 
determined by fit to some of the observables of the vacuum QCD. The various EOSs as obtained in 
the present work can be completely specified by (i) the source for the nuclear matter EOS as listed 
in Table IJ and (ii) the values of Go and Gy used in computing the EOS of the CSQM within the 
NJL model. The properties of spherically symmetric non-rotating and axially symmetric rotating 
compact stars are obtained by solving the Einstein's equations in ID and 2D, respectively. The 
numerical computations are performed by using RNS code written by Stergioulas and Friedman 
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The equilibrium sequence of compact stars for a given EOS is obtained by varying the central 
energy density e^- For the stable configuration, 

§^>0, (13) 

dec 

where, M is the gravitational mass of the non-rotating compact star. The equilibrium sequences 
for the non-rotating compact stars resulting from our EOSs are plotted as mass versus radius in 



Figs. |4]and[5l The central energy density increases as we move along these curves from the right 
hand side. The portion of the curves left to the solid circles represent the equilibrium sequences of 
hybrid stars with CFL quark matter core. The curves between the solid circle and triangle represent 
the sequences of hybrid star composed of 2SC quark matter core. It is clear from the lower panel 
of Fig. Hland upper panel of Fig. [5] that the stable configurations of hybrid stars with CFL quark 
matter core belong to third family of compact stars. Further, irrespective of the choice of the EOS 
of the nuclear matter, the stable configurations of the non-rotating hybrid stars exist within the 
NJL model only when the EOSs for the CSQM are constructed for Go = 1.2Gs with Gy < O.IG5. 
These values of Go and Gy, for which the stable configurations of the hybrid star exist, are very 
much similar to the ones found in Ref. jsS]. It appears that the stability of the hybrid stars with 
CSQM core depends solely on the choice of the EOS for the CSQM. However, the composition 
of the hybrid stars depend on the behaviour of the nuclear matter EOS. For instance, in case of 
the TMl and the NL3 EOSs of the nuclear matter, the core of the hybrid stars are composed of 
CSQM which is either in 2SC phase or in the CFL phase. In the later case, CFL quark matter core 
is surrounded by a layer of 2SC quark matter with the outer layer composed of nuclear matter . 
The thickness of the 2SC quark matter at the maximum hybrid star mass is around 0.5 -0.7km and 
its mass is ~ O.IMq. On the other hand, no 2SC quark matter appears in the stable configurations 
of the hybrid star constructed using the EOSs for which the nuclear matter part correspond to the 
APR and SLY4. 

The equilibrium sequences for the hybrid stars rotating with fixed rotation frequency / are 
constructed. As an illustration, in Fig. [6l we plot mass verses circumferential equatorial radius 
Req at fixed values of the rotational frequency obtained for two different EOSs. For the clarity, we 
mainly focus on the regions of the M - R^q curves corresponding to the sequences of the hybrid 
stars which are relevant in the present context. We see that beyond certain frequency, so-called the 
critical rotation frequency /cnt, the stable configuration for the rotating hybrid star does not exist. 
The solid black lines in Fig. [6]represent the result obtained at the / = /cnt. In Fig. |7]we plot the 
values for the /cnt (left panel) calculated for the cases for which the stable configurations for the 
non-rotating hybrid star exist. It is evident that the values of fcnt are quite sensitive to the choice 
of the EOS for the nuclear matter as well as the CSQM. Depending on the EOSs, /„,, varies in 
the range of 350 - 1275 Hz. We also plot in Fig. |7] (right panel) the maximum mass M^ax for 
the non-rotating hybrid stars with CSQM core. It is interesting to note from this Fig. that the 
values of Go and Gy for a given nuclear matter EOS can be so adjusted that the resulting hybrid 
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star has (a) the maximum mass in the non-rotating limit larger than 1.44Mo which is the most 
accurately measured value for the maximum mass of a compact star and (b) the maximum 
allowed rotation frequency is larger than the current observational limit of 716 Hz ^44]. 

In Fig. [8] we show the correlations between the values of the fcnt for the hybrid stars and the 
radius i?i 4 for the neutron star with the canonical mass. We see that /cnt is large if value of R14 
is also large. Thus, hybrid star constructed for a given EOS for the CSQM can rotate faster if the 
EOS for the nuclear matter is stiffer. The existence of the correlations between the values of /cnt 
and i?i 4 may be due to the fact that the pressure at which the nuclear to the quark matter transition 
occurs is closer to the values of P1.4 as can be seen from the lower and upper left panels of Figs. |2] 
and[3l respectively. 

Finally, we would like to compare the present results with coresponding ones obtained within 
the MIT bag model BTl kSn. The present results as obtained within the NIL model are significantly 
different with those for the MIT bag model. Within MIT bag model the EOS for the CSQM 
can be obtained by adjusting the value of the CFL gap parameter and the bag constant such that 
the resulting hybrid stars with CFL quark matter core are gravitationally stable upto the masses 
~ 2Mq in the static limit and the the maximum allowed rotation frequency is much larger than 
1 kHz. However, it can be seen from Fig. Ul stable configurations of the hybrid stars with CFL 
quark matter core obtained within the NIL model are having the maximum values for the mass 
and the rotational frequency appreciably lower than those obtained for the MIT bag model. The 
differences between the results for the MIT bag model and the NIL model can be attributed to the 
fact that the constituent quark masses, chiral condensates and the colour superconducting gaps in 
the later case are computed self-consistently as a function of baryon density. 



IV. CONCLUSIONS 



The stability of non-rotating and rotating hybrid stars, composed of the colour superconducting 
quark matter core surrounded by a nuclear mantle, is studied by using several EOSs. The EOSs 
for the nuclear matter, employed at lower densities, are based on the variational and the mean 
field approaches. We use a diverse set of nuclear matter EOSs such that the resulting maximum 
neutron star mass lie in the range of 2.2 - 2.8Mo and the radius at the canonical neutron star 
mass vary between 11.3 - 14.8km. The EOSs at higher densities corresponding to the colour 
superconducting quark matter in the 2SC or the CFL phase, are calculated within the NIL model 
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using different values of coupling strengths for the scalar diquark and isoscaler vector terms. The 
EOS at intermediate densities are obtained using a Maxwell construction. 

We find that the stability of the non-rotating hybrid stars is very much sensitive to the EOSs 
for the colour superconducting quark matter and almost independent of the choice for the EOS for 
the nuclear matter. The stable configurations of the hybrid stars exist only for the large enough 
value for the scalar diquark coupling strength. Though, the stability of the hybrid stars are not 
sensitive to the choice of the EOSs for the nuclear matter, but, compositions of the hybrid stars 
are at variance for these EOSs. If the EOS for the nuclear matter is stiff, core of the hybrid star is 
composed of colour superconducting quark matter which is either in the 2SC or the CFL phase. 
In the later case, CFL quark matter core is surrounded by a thin layer of the 2SC quark matter and 
the outer layer composed of nuclear matter . 

The stability of the rotating hybrid star is sensitive to the choice of the EOS for the nuclear 
matter as well as that for the colour superconducting quark matter. In particular, we find that 
the values of the critical rotation frequency vary from about 350 Hz to 1275 Hz depending upon 
the choice of the EOSs for the nuclear matter and the colour superconducting quark matter. Our 
results also indicate that the EOSs for the colour superconducting quark matter obtained within 
the NJL model may be adjusted for the various nuclear matter EOSs in such a way that it yields 
(a) the maximum mass in the non-rotating limit larger than 1 A4Mq which is the most accurately 
measured value for the maximum mass of a compact star and (b) the maximum allowed rotation 
frequency is larger than the current observational limit of 716 Hz. 

Finally, we would like to mention that our present study can be extended in several ways. The 
quark matter in the crystalline color superconducting phase, expected to appear at the intermediate 
densities, should also be considered. One might also include the contributions from the hyperons 
which would soften the hadronic EOS. The phase transition from hadron to the quark matter should 
proceed via mixed phase which can be constructed using the Gibbs conditions. 
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TABLE I: Values of the maximum mass Mmax and corresponding central energy density emax and radius 
Rmax obtained for different EOSs of the nuclear matter. The radius /?i 4 for the neutron star with canonical 
mass ( 1 AMq) are also given. 



EOS 


Approach 








R[ 4 


Ref. 






(10i5g/cm^) 


(Mo) 


(km) 


(km) 




APR 


Variational 


2.80 


2.19 


9.9 


11.3 


m 


SLY4 


NRMF 


2.84 


2.05 


10.0 


11.7 


[£7J 


BSRIO" 


RMF 


2.14 


1.97 


11.6 


13.3 


m 


TMl 


RMF 


1.87 


2.19 


12.4 


14.4 


m 


NL3 


RMF 


1.55 


2.79 


13.3 


14.7 





"This EOS is obtained using one of the several parameter sets of the extended RMF model given in our earlier work 
1 4811 ■ Each of these parameterizations corresponds to different values of the strength ^ for the a> -meson self-coupling 
term and neutron-skin thickness Ar in ^°^Pb nucleus. The remaining parameters of the models were calibrated to yield 
reasonable fit to the bulk nuclear observables and nuclear matter incompressibility coefficient. In the present work we 
use the parameter set with ^ = 0.03 and Ar - 0.2fm which will be referred henceforth as BSRIO. 
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FIG. 1: (Color online) The nuclear matter EOSs plotted as pressure versus energy density. The energy 
density is normalized by 60 = 150 MeV/fm^ which is the typical value of the energy density for the nuclear 
matter at the saturation density. 
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FIG. 2: (Color online) The pressure as a function of the quark chemical potential (left panel) and the baryon 
density (right panel) for the nuclear and quark matter. The EOSs for the quark matter in the 2SC and CFL 
phases are obtained within the NJL model using Gy = with Go = I IG5 and I.2G5. The solid circles 
on the various EOSs for the nuclear matter indicate the pressure at the center of the neutron star with the 
canonical mass (IAMq). 
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FIG. 3: (Color onUne) Same as Fig. El but, for Gy = O.IG5 and O.2G5 with Go - 1.2Gs. 
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FIG. 4: (Color online) Plots for the mass-radius relationships for the equilibrium sequences of non-rotating 
compact stars obtained using various EOSs as shown in Fig. |2l The curves on the left of the solid circles 
represent the equilibrium sequences for the hybrid stars with core composed of the quark matter in the CFL 
phase. The curves between the solid circles and triangles represent the hybrid stars with 2SC quark matter 
core. The absence of solid triangle on a curve means that the hybrid stars contain quark matter only in the 
CFL phase. 
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FIG. 5: (Color online) Plots for the mass-radius relationships for the equilibrium sequences of the non- 
rotating compact stars obtained using various EOSs as shown in Fig. [3l The solid circles and triangles 
divide the curves according to the composition of the hybrid stars as described in Fig. |4l 
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FIG. 6: (Color online) Plots for the mass verses circumferential equatorial radius /?eq at fixed values of the 
rotational frequency as indicated along each of the curves (in Hz). The black solid lines represent the results 
obtained at the critical frequencies faw- For / > /crit, the stable configurations of hybrid star do not exist. 
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FIG. 7: (Color online) The values of the critical rotation frequency fan (left panel) for the hybrid stars with 
the CSQM core and their maximum masses (right panel) in the non rotating limit as obtained using different 
EOSs for the nuclear matter and the CSQM. 
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FIG. 8: (Color online) Correlations between values of /oit for the hybrid stars and the radius R14 for the 
neutron star with canonical mass as listed in Table Jl 
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